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We study the effects that the gravitational interaction of f(R) theories of gravity in Palatini 
formalism has on the stationary states of the Hydrogen atom. We show that the role of gravity 
in this system is very important for lagrangians f(R) with terms that grow at low curvatures, 
which have been proposed to explain the accelerated expansion rate of the universe. We find that 
new gravitationally induced terms in the atomic Hamiltonian generate a strong backreaction that 
is incompatible with the very existence of bound states. In fact, in the 1/R model, Hydrogen 
disintegrates in less than two hours. The universe that we observe is, therefore, incompatible with 
that kind of gravitational interaction. Lagrangians with high curvature corrections do not lead to 
such instabilities. 
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I. INTRODUCTION 

The accelerated expansion rate of the universe [l| is 
one of the biggest puzzles that theoretical physics faces 
nowadays. Dark energy sources within the framework 
of General Relativity (GR) have been postulated as the 
missing element that could explain that phenomenon. 
On the other hand, modified theories of gravity have 
been proposed as an alternative to dark energy sources. 
Modified theories usually provide self-accelerated cosmic 
solutions on purely geometrical grounds, making unnec- 
essary the introduction of unobserved exotic sources of 
matter-energy. The modification of the gravitational 
laws is, however, a very delicate issue in which intuition 
is not always a good guide. In fact, a modification 
originally thought to affect the dynamics at large 
scales could end up having non-trivial effects at shorter 
scales. It is thus necessary to study the dynamics of 
the different gravitational theories in different regimes 
and identify their positive and negative aspects aiming 
at learning how to construct theories that exhibit the 
desired properties. In this sense, theories of gravity in 
which the lagrangian is some function f(R) of the scalar 
curvature R manifest many interesting properties and 
have attracted much attention in the recent literature. 

The equations of motion of f(R) theories can be 
derived in two different ways depending on whether the 
connection is seen as independent of the metric (Palatini 
formalism) or as dependent of it (metric formalism). In 
the metric formalism, i.e., when the connection is the 
Levi-Civita connection of the metric, besides the metric 
one identifies an additional scalar degree of freedom, 
which turns the scalar curvature R into a dynamical 
object. The interaction range of this scalar field depends 
on the form of the lagrangian and can change due 
to different reasons. When perturbation theory is 
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applicable, the first order approximation shows that, for 
models of interest in the late-time cosmic acceleration, 
the interaction range changes driven by the cosmic 
expansion The scalar field can be short-ranged for 
some time (radiation and matter dominated eras) and 
then turn into a long-ranged field, causing late-time 
cosmic acceleration. This type of lagrangian is ruled out 
by solar system experiments. Some f{R) lagrangians, 
however, cannot be treated perturbatively. Nonetheless, 
one can still define an effective mass or interaction 
range for the scalar, which now depends on the local 
matter density. The constraints on such lagrangians by 
local experiments have been discussed recently in the 
literature 3. 

In this paper we will consider the other formulation of 
f(R) theories, namely, the Palatini formalism. Almost 
surprisingly, allowing the connection to be determined 
by the equations of motion does not introduce new 
dynamical degrees of freedom. In these theories, the 
metric turns out to be the only dynamical field, which 
satisfies second-order differential equations. As we will 
see, the effect of the lagrangian f(R) is to change the 
way matter generates the metric by introducing on the 
right hand side of the field equations new matter terms 
that depend on the trace of the energy-momentum 
tensor of the sources. In vacuum, the field equations 
reduce (always and exactly) to those of GR with a 
cosmological constant. For this reason, it has been 
thought for some time that Palatini f(R) theories could 
pass the solar system observational tests 0, Q (see 
also [6[ for a discussion of this point). It has also been 
shown that different choices of lagrangian f(R) are able 
to accommodate several of the different cosmic eras of 
the standard cosmological model [7]. However, a more 
careful analysis of the gravitational dynamics in the 
presence of sources indicates that these theories might 
be in strong conflict with our understanding of the 
microscopic [23[ world [1, H, [l(| (the study of polytropic 
matter configurations also points in this direction pl|). 
This aspect, together with the cosmological viability 
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of different models, is of primary importance due to 
the fact that almost all cosmological observations relay 
on the detection of electromagnetic radiation, which is 
intimately related to the quantum mechanical nature 
of atomic and molecular structure. In this work, we 
elaborate in this direction and study how the gravita- 
tional interaction in Palatini f(R) theories affects the 
non-relativistic limit of the (one-particle) Dirac equa- 
tion. The analysis of various gravitationally-induced 
correcting terms in the resulting Schrodinger-Pauli 
equation will provide us with solid arguments against 
the existence in the gravity lagrangian of correcting 
terms relevant at low cosmic curvatures. As we will see, 
the very low matter density (or curvature) scales that 
characterize (infrared-corrected) modified lagrangians 
can be reached near the zeros of the wavefunctions. This 
causes a strong gravitational backreaction that makes 
unstable the stationary states of the Hydrogen atom. In 
particular, we find that the ground state disintegrates 
in a matter of hours. On the contrary, if the gravity 
lagrangian is modified by high curvature corrections, the 
backreaction effects are negligible and the atom remains 
stable. 

The paper is organized as follows. In section II, we de- 
fine the action of f(R) theories in Palatini, derive the field 
equations, and discuss the metric generated by micro- 
scopic systems. We then introduce two illustrative mod- 
els, namely, f(R) = R + -jf- and f(R) = R - which 
will help us compare the fjehavior of the metric when 
the GR action gets ultraviolet or infrared corrections, re- 
spectively. In section III we derive the non-relativistic 
limit of Dirac's equation starting from its curved space- 
time formulation. We then discuss the effects induced by 
the modified Schrodinger-Pauli equation in the station- 
ary solutions of the Hydrogen atom. We conclude with a 
summary and discussion of the results obtained. In the 
Appendix we estimate the decay rate of the ground state 
of the atom. 



II. THE THEORY 

Let us begin by defining the action of Palatini theories 

%, T, Vm] = ^2 J <^ x V-9f( R ) + Sm [g»u, 1pm] (1) 

Here f(R) is a function of R = g^R^T), with Rn V (T) 
given by R^(T) = -d^l + d x T^ + I*,r' A - T^T^ 

where r^„ is the connection. The matter action S m de- 
pends on the matter fields ip m , the metric g^, which 
defines the line element ds 2 = g fiu dx^dx u , and its first 
derivatives (Christoffel symbols). The matter action does 
not depend on the connection T*, which is seen as an 
independent field appearing only in the gravitational ac- 
tion (this condition is not essential and can be relaxed at 
the cost of introducing a non- vanishing torsion) . Varying 



(fTj) with respect to the metric <7 M „ we obtain 



f'(R)R pi/ (r) - -f{R)g, v = n 2 T^ u 



(2) 



where f'(R) = df /dR. From this equation we see that 
the scalar R can be solved as an algebraic function of the 
trace T. This follows from the trace of (J2j) 



f'(R)R-2f(R) = K, 2 T, 



(3) 



The solution to this algebraic equation will be denoted 
by R = K(T). The variation of © with respect to T* v 
must vanish independently of @ and gives 



= (4) 



where /' = f'(1Z[T]) is also a function of the matter 
terms. This equation leads to 



d v tp 



(5) 



where t. 



and 



"\w = Ybtiu, v = /'(TC[n]j is dimensionless 
and normalized to unity outside of the sources (T = 0) . 
It is now useful to rewrite adding and subtracting 

^K(T)g^ = ^R a p{T)t^ to get 



[Kf'- f\. 

2(f> 



(6) 



where G M „(i) is the Einstein tensor associated to t^ v . 
The equations of motion ^ for the auxiliary metric 
are considerably simpler than those for g^ Vl 



R^(g) - -^g^ v R{g) 



T 



Kf'- f 
2/' 



f 



2(/') 2 

+ jl^^uf -g^nf] , (7) 

because of the difficulty of dealing with the matter 
derivatives df ~ dT and d 2 f ~ (dT) 2 ~ (d 2 T). 
Solving for using the system ([6]) and then go- 
ing back to g^ via the conformal transformation 
9i±v = 4>{T)~ 1 t^ l , is a useful simplification that makes 
the task of finding solutions much easier. Note that 
this fortunate circumstance is due to the fact that the 
conformal transformation completely cancels out the 
disturbing derivatives of 0. Moreover, the conformal 
relation between the two metrics puts forward the fact 
that the metric <? M „ receives two kinds of contributions: 
non-local contributions that result from integration 
over the sources, which produce the term t llVl and 
local contributions due to <f>(T), which depend on the 
local details of T at each space-time point. This local 
contribution arises due to the independent character 
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of the connection and, to our knowledge, does not 
appear in any other metric theory of gravity, where the 
connection is generally assumed to be metric compatible 
(Levi-Civita connection). 

The similarity between the field equations of GR and 
© suggests that for weak sources and reasonable choices 
of lagrangian f(R) (those that lead to negligible cos- 
mological term )» the right hand side of ([B]) is 

small and, like in GR, t^ v can be expressed as ta V {x) ~ 
tlnv + h IJlU (x), where ^^(x)! <C 1 is given as an integral 
over the sources (for details of exact calculations see [8j]). 
For very weak sources such as atoms or elementary par- 
ticles, we find that the self-contribution to — > 0. On 
the other hand, if our microscopic system is placed in an 
external gravitational field whose contribution to t^ is 
not negligible, we can always take a coordinate system in 
which the metric at the boundaries of a box containing 
the system (large as compared to the microscopic system 
but small as compared to the range of variation of the 
external metric t^ u ) becomes ~ ?7 M „. In both situations, 
the metric g M „ becomes simply 

g»u(x) « <KT)~ V (8) 

where <f>(T) — > 1 in the boundaries of our auxiliary 
box but might depart from unity within the sources, 
depending of the particular lagrangian chosen. 
If one trivializes the role of the local term 4>(T), then 
one finds that relative motion between particles is not 
very much affected by the modified gravity lagrangian 
(T3 |. However, we find that the presence of (/((T)" 1 in 
front of r/fiu is very important because the matter fields 
in (fTJ) are coupled to g^ and, unlike in all other known 
metric theories of gravity, g^ v only becomes locally 
•q^v in regions where T = exactly. Consequently, the 
4>(T) dependence of g^ v induces new interactions and 
self-interactions between the matter fields [8j, as will 
be explained here in detail. Note that the presence of 
the (scalar) term <f>(T) is physical and not a problem 
of choosing the wrong coordinate system, as criticized 
in [l3|. In fact, if one computes geometrical invariants 
such as 24J R a j3 1 x R a ^ 1 \, various derivatives of <f>(T) 



appear and cannot be eliminated by choosing different 
coordinate systems because R a j3 1 x R a P 1 \ is a coordinate 
invariant. Note also that the dependence of R a py X R a P J \ 
on the local energy-momentum distribution via <fi(T) 
tells us that the geometry might be subjected to mi- 
croscopic fluctuations driven by the fluctuations of T 
and modulated by the form of the gravity lagrangian 
[recall that 0(T) = f'(K[T])/f'(K[0])]. Therefore, 
lagrangians sensitive to low energy scales could lead to 
unnacceptable microscopic curvature fluctuations, while 
others could lead to more robust geometries which would 
only fluctuate at very high energies. In this latter case, 
however, neglecting the contribution of h^ v could not 
be well justified (think for instance in the hypothetical 
production of black holes in particle accelerators) . 



A. Two illustrative f(R) models. 

We will now study the behavior of the function (j>(T) 
f'[TZ(T)]/f'[TZ(0)] for two illustrative models. 



1. Ultraviolet corrections: f(R) = R + 

This model is characterized by a high energy/curvature 
correction R 2 /Rp, where the subscript P stands for 
Planck scale. In this case, we find that 1Z(T) — —k 2 T is 
the same as in GR, and the function 4>(T) is given by 



<f>{T) = 1-2 



K 2 T 

Rp 



(9) 



We thus see that only at very high matter/energy den- 
sities will the function 4>(T) significatively depart from 
unity. 



2. Infrared corrections: f(R) 



This model was initially proposed in [l4( within the 
metric formalism and is characterized by a low curvature 
scale fx 2 . In this case, we find that 1Z(T) = — (k 2 T + 
yJ{n 2 T) 2 + 12^ 4 ) /2 recovers the GR limit for |k 2 T| > p 2 
and tends to a constant 1Z ~ p 2 for |k 2 T| <C p 2 ■ The 
function </>(T) is given by 



i 



i 



2[1 + + 12/r 2 



(10) 



Here r = -T/T c , T c = /i 2 /k 2 = Pfi , and 
Pfj, ~ 10 -26 g/cm 3 represents the characteristic 
cosmic density scale of the theory, which triggers the 
cosmic speedup. It is easy to see that at high densities, 
as compared to p M , <f>(T) — > 3/4, whereas for p <C 
we find <p(T) — > 1. Note that we could have chosen the 
normalization of <f>(T) differently and in such a way that 
at high densities 4>(T) — > 1 whereas at low densities 
4>(T) — * 4/3. In this latter case, however, the physical 
metric g^ v would tend to §?7 M „ in vacuum. We find 
our first choice a more natural normalization (though 
arbitrary anyway), since it makes g^ v = rj^ in vacuum. 



III. DIRAC EQUATION IN CURVED SPACE 

It is well known 

ee nam 

that the energy levels of a 
Hydrogen atom falling freely in an external gravitational 
field (in GR) will be shifted in a very characteristic way 
due to the interaction of the electron with the curvature 
of the space-time. Though external fields in Palatini the- 
ories of gravity must also lead to this phenomenon, we 
will focus here on a different aspect. We will study the 
effect that the local energy-momentum densities have on 
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the non-relativistic limit of the Dirac equation due to the 
factor ^(r)" 1 appearing in (jHJ). 

For the sake of clarity, let us briefly consider the different 
contributions that make up T, which generate the met- 
ric ([8]) seen by the system. The electromagnetic field, 
which is treated as classical, is traceless and, therefore, 
does not contribute to T. The atomic nucleus can be 
modeled as point-like or as described by an extremely lo- 
calized wave-packet contributing with T/v = —mjsrS e (x), 
where S e is some representation of the Dirac delta func- 
tion with spread e centered at the origin, and mpj is the 
nuclear mass. The motion of the electron is described 
by the one-particle Dirac equation, which generalized to 
curved space-time [HI, [H, [13] can be derived from the 
following action (the notation will be explained below) 



S m [g^,Tp] 



(11) 



Upon variation of this action with respect to g^ v one finds 
the energy-momentum tensor associated to the electron, 
whose trace is given by [l8| 

T e = -m^ip (12) 

In summary, T = Tn + T e = —rriNS e (x) — rmpip. 



A. Derivation of the non-relativistic limit 

From the action (TTTj) , we can derive the curved space- 
time version of Dirac's equation 



{iX^Df, - m)V> = 



(13) 



Here A p = e%j a are the curved space Dirac matri- 
ces, which are related to the constant Dirac matrices 
{7a, 7b} = 2i]ab by the vierbein (recall that g^" = 
W e a e b)' The covariant derivative is given by 



Even though ([15]) is not, in general, completely separable 
due to the non-linearities introduced by the dependence 
of T on ipip, stationary solutions do exist. To find them, 
it is useful to write the equation in the form idtip = Hip 
[25j as follows 

id t ip= a ■ (p- eA + iVfi) + (eA + id t Q) + m(3 tp 

Let us now focus on the positive energy solutions of this 
equation. It is easy to see that taking ip(t, x) — e~ iEt £(x) 
we have T e = — m££, dt^l = 0, and (TTg|) turns into 



E£ = [a ■ 7? + eA a + fh(3\ £ 



(19) 



where we have used the shorthand notation ir = (p — 
eA + iVQ), Denoting by r\ and \ l ne large and small 
components, respectively, of the Dirac spinor 



we find the following relations 
1 



X = 



E- 



Et) 



eA 
1 



a ■ TTT] 



a ■ 7T- 



E- 



eA 



eA 



[to + E - eA)] 2 



(20) 

(21) 
V (22) 
77(23) 



We will now proceed to compute the lowest-order non- 
relativistic limit. We first decompose the energy E in two 
parts, E — mo + £, where too is a constant of order ~ m 
(to be discussed further below) and £ <C 777 represents 
the non-relativistic energy. We then expand assuming 
that the rest mass is much larger than the kinetic and 
electrostatic energies, to ~ too 3> \£ — eAo\, and retain 
terms only of order 1/toq. The above relations reduce to 



D„ 



d u + ieA u 



(14) 



with w^ b representing the spin connection, A^ is the 

electromagnetic vector potential, and Y, a b = j[7<u7&]- 
Since, by construction, the matter action is not cou- 
pled to the connection r" , the spin connection w® b 
must be defined in terms of the Christoffcl symbols 

c uu = \- {®v,9pv + dvdpv - dp9n») and tnc vierbein as 



w 



nb 



e^V tl e ab = e a a {d^e ab + C° x e xb ). From © it is 

easy to see that e° = (p^ 1 ! 2 ^ and e% = ^^S^. After a 
bit of algebra, (TT3)) turns into 



[i^ida + ieA a - d a n) - m] V = 



where we have defined 



fl = (3/4) ]n<j)(T) 



(15) 

(16) 
(17) 



X 
£y 



1 



- too 
1 



-a- ttt) (24) 



too 



(a ■ tt) 2 + (to - toq) + eA 77 (25) 



-rrnf [I- 0(|t?| 2 /too)] 77 = -mrfrj (26) 



The wavefunction of the electron is then identified with 
77, which to this order coincides with the positive energy 
Foldy-Wouthuysen bispinor 19]. From (j2l))) we see that 
the non-linearities contained in 7? and 777 in Q25|) only 
depend on 77. Expanding the operator (a ■ if) 2 we find 



£ V = 



1 



to + too 
1 



-[(p- eA) 2 - eS ■ B] + eA \ 77 



x V) - 2ie(A ■ VO) (27) 

v 2 n- \vn\ 2 + 2(vn- v) + (to-to )}?7 



to + to 

2n IV70I2 
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The first line of this equation is very similar to the 
well-known non-relativistic Schrodinger-Pauli equa- 
tion (see (f29| below). The only difference being 
the term 1/(to + mo). The second and third lines, 
however, represent completely new terms generated 
by the Palatini gravitational interaction. When the 
gravity lagrangian is that of GR, 4>(T) = 1, we recover 
the Schrodinger-Pauli equation if toq is identified with m. 



IV. APPLICATION: THE HYDROGEN ATOM 

To gain some insight on the role and properties of the 
various terms in (|2~T|) . we will proceed as follows. We first 
solve (EH) in the case of GR, f(R) = R, cf>(T) = 1, which 
is well known. Then we switch to a different gravity 
lagrangian (assuming that we have the ability to do 
that) and study how the system reacts to that change. 
The reason for this is that in a general f(R) the metric is 
sensitive to the local T M „ via ^>(T) -1 , and changes in the 
metric due to the matter distribution could react back 
on the matter equations. If the new interaction terms in 
(|2T|) lead to small perturbations, then the initial wave- 
functions will be, roughly speaking, stable with perhaps 
small corrections which could be computed using stan- 
dard approximation methods. If, on the contrary, the 
energy associated to the gravitationally-induced terms is 
large, that would mean that the original configuration is 
not minimizing the modified Hamiltonian and, therefore, 
large modifications would be necessary to reach a new 
equilibrium configuration. Depending on the magnitude 
of the reaction on the system, we could estimate whether 
the theory is ruled out or not. 

Let us first consider the f(R) model with ultraviolet 
corrections introduced in section Hi A II In this case, the 
function 4>{T) = 1 — 2-^-^ can be expressed as 



2[p N (x) + p e (x)} 



Pp 



(28) 



where pn(x) = 
P e (x) = rf(x)r](x 
p P = Rp/k 2 
(Planck scale). 



m N 5 t {x), pe(x) = mP e (x), 
is the probability density, and 
is a very high matter-density scale 
Since the scale pp is much larger 
scale 



than any density scale reachable by the electron 
wavefunction and even by the very peaked nu- 
clear wavefunctions (pn/pp ~ 10 -79 ), we see that 

f2(T) « 3 PN(x)+ Pc (x) ^ ^ m(1 _ PN ( X )+ Pc (x) j 

lead to strongly suppressed contributions (in fact, they 
are much smaller than the corresponding Newtonian 
corrections \hoo\ = GM/c 2 Rn ~ 10~ 39 ). Identifying too 
with the electron mass to, the leading order corrections 
to the wave functions and the energy levels could be 
computed by perturbation methods and would lead to 
virtually unobservable effects. 



Let us now focus on the model with infrared correc- 
tions introduced in section lll A 21 Expressing length units 
in terms of the Bohr radius (eto ~ 0.53 • lCP 10 m), we 
find t = pe ^ = 10 24 P P (a;), where we have intention- 

ally omitted the nuclear contribution (only relevant at 
the origin) for simplicity. This expression for r indicates 
that the electron reaches the characteristic cosmic den- 
sity, t ~ 1, in regions where the probability density is 
near P e {x) ~ 10~ 24 . In ordinary applications, one would 
say that the chance to find an electron in such regions is 
negligible, that that region is empty. In our case, how- 
ever, that scale defines the transition between the high 
density (t ^S> 1) and the low density (t <C 1) regions. 
In regions of high density, we find that (j) rapidly tends 
to a constant, = 3/4, which leads to to = 2to/\/3 
and Vfl = 0. If we then identify to — > a/3too/2, equation 
(|2T|) reduces to the usual Schrodinger-Pauli equation 



I 2m 



[(p - eAf - ea ■ B] + eAo \ r\ 



(29) 



This fact justifies the introduction of too above. Let us 
now see what happens in regions of low density. In those 
regions, <fi(T) tends to unity, Vf2 = 0, and to — > to as 
t — * 0. The mass factor dividing the kinetic term is 
now a bit smaller (too > to) than in the high density 
region. But the mass difference m — too is not zero. This 
is a remarkable point, because to — too ~ — 0.13too is 
negative and of order ~ too, which represents a large 
contribution to the Hamiltonian. To better understand 
the effect of this term, it is useful to consider the ground 
state, 



77(1,0,0) 



where | i , s) represents a 



normalized constant bispinor. In this case, the transition 
from the high density region to the low density region 
occurs at r w 26cto. In FigfTjwe have plotted the most 
representative potentials in dimcnsionless form 



Vn 



2tooC o,q 



to 



-d x n + din - \d x ii\ 

x 



(30) 
(31) 
(32) 



where V e is the electrostatic potential generated by the 
proton, lengths are measured in units of the Bohr ra- 
dius, x = r/ao, and energies in units of 2 m a' 1 ~ ^3-6 

eV. Note that Vn = V 2 f2— |Vf2| 2 only contains the most 
important contributions associated to SI. In this case, 
Vh represents a small transient perturbation. The mass 
difference V m , however, introduces a deep potential well 
in the outermost parts of the atom that must have im- 
portant consequences for its stability. (Note that this 
effect is not an artifact of the non-relativistic approxi- 
mation, since it also occurs in the full relativistic theory 
(IT5|) - (|T5]) due to the density dependence of to). In the 
initial configuration of the atom, corresponding to GR, 
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Figure 1: Contribution of the different potentials in the 
ground state. The solid line, which represents the sum of 
all the potentials, tends to the constant value — 0.134mo (or 
—5048 in the units of the plot). 



Figure 2: The different contributions in this plot are V e ~ —1, 

V m 5 • 10 3 , and Vh ~ ±10 19 . The y-axis is measured in 

units of 13.6eV; the x-axis in units of ao. 



the wavefunction of the ground state is concentrated near 
the origin, where the attractive electric potential is more 
powerful (V e — > — oo). As we switch on the 1/R theory, a 
deep potential well of magnitude ~ — 0.13m appears in 
the outer regions of the atom, where p e {x) < p M , which 
makes the ground state unstable and triggers a flux of 
probability density (via quantum tunneling) to those re- 
gions. The half life of Hydrogen subject to this potential 
can be estimated using time dependent perturbation the- 
ory (see the Appendix) yielding 

t = - w 6 • 10 3 s (33) 

We thus see that the initial, stable configuration is 
destroyed in a lapse of time much shorter than the age of 
the Universe, which is in clear conflict with experiments 

Further evidence supporting the instability of the 
atom is found in the existence of zeros in the atomic 
wavefunctions in between regions of high density be- 
cause, obviously, before (and after) reaching p e (x) = 
the characteristic scale p e (x) ~ p^, is crossed. The first 



excited state, 77(2,0,0) = 




has a zero at r = 2a$. The radial derivatives of 4>(T) 
at that point are very large and lead to very important 
perturbations which overwhelmingly dominate over 
any other contribution (see Fig. 2). The magnitude of 
Vn = V 2 ft - |Vf)| 2 at r = 2a oscillates between 10 20 
and — 10 21 eV in an interval of only 2 • 10 -10 ao. Needless 
to say that this configuration cannot be stable and that 
strong changes must take place in the wave function to 
reduce the energy of the system. Such changes should 
tend to reduce the magnitude of the density gradients 
(Vf2) to minimize the value of Vq, which will likely lead 
to a rapid transition to the ground state, where Vn is 
small. One can easily verify that strong gradients Vfi 
also appear at the zeros of all the ?7n,o,o wavefunctions, 



which generate large contributions Vh in those regions. 
Furthermore, if one considers stationary states with 
I ^ 0, Vq has important contributions not only at the 
zeros of the radial functions, but also at the zeros of 
the angular terms. Thus, the pathological behavior 
described for the spherically symmetric modes gets 
worse for the I ^ states. One thus expects the decay 
of these states into states with less structure (weaker 
gradients) such as the ground state, which will later 
decay into the continuum. All this indicates that the 
existence of bound states, with localized regions of high 
probability density (where "high" means above the scale 
Pij), are impossible in this theory because of the large 
gradient contributions Vh and the deep potential well V m . 



V. SUMMARY AND DISCUSSION 

In this work we have deepened into the effects that 
the matter-energy density dependence of the metric in 
Palatini f(R) theories has at microscopic scales. In 
particular, we have studied the effects that switching 
from GR to a different gravitational interaction, such as 
the f(R) = R + R 2 /R P or f(R) = R-p 4 /R models, has 
on the stationary solutions of the Hydrogen atom. To do 
so, we started with the Dirac equation in curved space 
and computed its non-relativistic limit. Then we looked 
at the contribution of the different new interaction 
terms appearing in the resulting (effective) Hamiltonian 
(|27p . We have found that the existence of bound states 
in theories with infrared corrections is problematic 
for several reasons. Firstly, due to the dependence 
of the effective electron mass fh — m(j>~ 2 (T) on the 
matter-energy density T = —m^jt/j, the effective mass 
seen in the inner (high density) regions and the outer 
(low density) regions of the atom is not the same. This 
generates a potential well that triggers the tunneling of 
probability density from the inner parts to the outermost 
parts of the atom, which eventually disintegrates the 
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atom. Secondly, wilder perturbations arise in those 
points and directions in which the wavefunction has 
zeros. This is due to the contribution of terms like V 2 f2 
and |V51| 2 when the characteristic scale p^ is crossed 
[27! ]. Minimizing the contribution of those terms would 
require a transition to states with less pronounced 
gradients such as the ground state, which would latter 
disintegrate into the continuum via tunneling. 

Though these instabilities have been discussed within 
the non-relativistic limit, we do not find any reason to 
attribute their existence to an artifact of this approxi- 
mation. In fact, the dependence of the mass on the local 
energy density was already apparent in (fl"5|) - (fl"8|) . In 
addition, derivatives of <fi(T) appear in the term c^f^T). 
Therefore, the relativistic description seems unable to 
cure the pathologies found in the non-relativistic limit. 
In addition, one can also check, by direct calculation 
of R af}l x R a ^ x ~ (d 2 fl) 2 + . . ., that the space-time 
geometry is strongly fluctuating and far from being fiat 
in those regions where p e ~ p^. 

Our results are also likely to hold even in the case in 
which the spin connection in the matter action is kept 
independent of the metric. In that case, the connection 
has a non-vanishing torsion, though the metric (and 
hence the vierbein) is still conformally related to the 
metric associated to the connection (see |20j), which 
is the key to get terms of the form rrufr 1 / 2 and d^fl. 

Though we have only analyzed in detail the infrared- 
corrected model f(R) = R — p 4 /R, the instabilities asso- 
ciated to the potential well to — too and the zeros of the 
wavefunction must be present in all gravity models sen- 
sitive to low curvature/energy-density scales. Since the 
matter, as we know it, would be unstable in those theo- 
ries, the cosmological models considered in that context 
are empty of significance (see Q for a list of references). 
On the contrary, models which introduce deviations from 
GR at high curvatures, such as f(R) = R + R 2 /Rp, do 
not have any relevant effect on the atomic structure if 
the characteristic scale is sufficiently high. To reach and 
excite the high energy-density scale one should deal with 
highly localized wave-packets, which will surely require 
the consideration of quantum fields. The quantization of 
the matter fields then opens an exciting window to new 
phenomena. In fact, when -0 is seen as a quantum field, 
the function T appearing in ([8]) and (fl~5|) - (fl"8|) must be 
interpreted as (T), i.e., the quantum expectation value 
of the operator T in a given state. The Hamiltonian of 
the theory then depends on the particular quantum state 
under consideration through the expectation value (T). 
A direct consequence of this is that the time evolution of 
the states in the Hilbert space of the theory is nonlinear 
(2l| . This highly non-trivial fact could be used to impose 
tight constraints on the form of the gravity lagrangian in 
Palatini theories via quantum experiments. In fact, we 
believe that in order to guarantee the linear evolution of 
quantum states, it could be necessary that the gravity 



lagrangian were exactly that of Hilbert-Einstein. 
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Appendix 

We briefly sketch here the computation of the half 
life given in (|33|) . Our calculation will be approximate 
and should provide a reasonable estimation of the or- 
der of magnitude of r = h/T. We will first assume 

that the kinetic term 



-=-7 — V can be approximated 

by — everywhere, even though to + too rs 1.87too 

in the low-density regions (recall that to = to/ \J<j>(T) 
becomes too when \T\ — > 00). Secondly, we will neglect 
the contribution of Vh and will approximate V m (r) 

by a step function of magnitude W$ = 



to 
too 



0(°°) 
0(oc) 



0(0) 



s — 0.13too in the region r > 26ao and 

zero elsewhere (see FigfT]). The total potential (for / = 0) 
when the l/R interaction is turned on can thus be seen 
as 



V(r) 



if r < 26ao 
-0.13mo if r > 26ao 



(34) 



This way we have reduced our problem to that of an 
initially stable bound state that becomes unstable and 
decays into the continuum when the initial potential 
U(r) — — - is transformed into V(r). This simplified 
scenario captures the essential features of our problem. 
The decay rate can be estimated using time-dependent 
perturbation theory. A simple and compact expression 
for the width T of a quasistationary state (which initially 
was a true bound state) is given by the following formula 
(see [lU for details) 



Ak2 2 

— — \MR)xk(R)\ 2 



(35) 



In our case, a — 1/ao, h = -\/2too(0.13tooc 2 — |e|)/fi, 
e = — 13.6eT^, t^o(R) represents the radial part of the 
partial wave expansion of the ground state evaluated 
at R = 26a , MR) = J hf^~ R ' a \ Xk{r) represents 
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the outgoing continuum mode, and Xk(R) = i a ° k 2 ■ version of the 1/R theory would disintegrate in less than 

Putting these numbers in (|35p . we find (|33p . which im- 
plies that the ground state of Hydrogen in the Palatini 
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